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Abstract. - We use computer simulations to compare the dynamical behaviour of torus and 
even- twist knots in polymers under tension. The knots diffuse through a mechanism similar to 
reptation. Their friction coefficients grow linearly with average knot length for both knot types. 
For similar complexity, however, the torus knots diffuse faster than the even twist knots. The knot- 
length auto-correlation function exhibits a slow relaxation time that can be linked to a breathing 
mode. Its timescale depends on knot type, being typically longer for torus than for even-twist 
knots. These differences in dynamical behaviour are interpreted in terms of topological features 
of the knots. 



The scientific study of knots has a long history, dating 
back at least to Johann Friedrich Gauss in the early 19th 
century. More recently, deep connections between knot 
theory, statistical mechanics [1] and quantum field the- 
ory [2] have stimulated a great deal of research both in 
physics and mathematics. 

The discovery of knots in bacterial DNA [3] and pro- 
teins [IHS] , as well as earlier work on synthetic polymers 
helped broaden the relevance of knot theory to chemistry 
and biology. Indeed, knots and links can be introduced in 
vivo into cellular DNA through processes such as replica- 
tion [5] and are regulated through enzymes such as type II 
topoisomerases that can both knot and unknot DNA 0. 
■ Biologists have exploited these topological effects to make 
I many discoveries about the nature of DNA in cells [10] . 
DNA is thought to be highly knotted inside some viral 
capsids [n], and we have recently suggested that these 
knots may control the ejection speed of a bacteriophage's 
DNA [H]. 

The discovery of knots in nature raises further ques- 
tions about the dynamics of knotted polymers. In con- 
trast to unknotted polymers, where dynamic behaviour is 
fairly well understood [13j . many basic questions remain 
open. In polymers without tension, a long topological 
timescale, not predicted by the standard Rouse model pTSl, 
was first discovered in the relaxation of the radius of gy- 
ration for knotted polymers [13] . Later work showed that 
this topological relaxation time decreased with increasing 
knot complexity pK], and could also be observed by mea- 



suring the knot length autocorrelation function [16]. Re- 
lated work [TTlfTS] showed that, by contrast, this timescale 
increased with knot complexity when knotted ring poly- 
mers are cut. 

Knots can be artificially introduced into biopolymers 
such as actin filaments [19] and linear DNA [20] with sin- 
gle molecule techniques. Bao et. al. [50] measured the 
diffusion coefficients of DNA knots by fluorescence mi- 
croscopy. The knot friction coefficient grew roughly lin- 
early with knot length, so that more complex knots dif- 
fused slower. Similar diffusion coefficients were obtained 
with simulations by Vologodskii [21]. In these studies 
polymers were under tension. It has also been argued that 
knots could exist in metastable tight states due to entropy 
alone [22] . 

The experiments performed in ref. [20] utilized tensions 
0.1-2 pN, which may be in a similar range to forces in- 
duced in vivo by enzymes such as polymerases [23] . Huang 
and Makarov |24j showed that the experiments were in 
the intermediate "elastic" regime where the external force 
/ > /csT/Zp, with Ip the persistence length. Here the 
force aligns the segments of the chain in the general di- 
rection of the force, and the knot size is determined by 
the bending elasticity of the chain vs the force. The diffu- 
sion coefficient is only weakly dependent on tension, an ef- 
fect also observed in experiments [20]. The elastic regime 
is bounded from above by a "tight knot" regime where 
/ ^ ksT/lp and the knot properties are dominated by 
details of molecular interactions between monomers. At 
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much lower tensions / ^ ksT/lp in the "blob" regime, 
the knot size fluctuates widely. 

Most previous work on knot dynamics has focussed on 
the "blob" regime. However, experiments on artificially 
knotted DNA [5D], as well as the expectation that some 
knots in vivo may be held under tension [12II23| , mean it is 
important to explore the effect of topology on dynamics in 
the elastic regime. This task will be our main focus here. 
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Fig. 1: (a) The first three torus knots, 
even-twist knots. 



(b) The first three 



We focus on two classes of knot topologies: torus knots 
and even- twist knots. The first few knots in each group 
are shown in Fig.[T]- knots of higher complexity have more 
crossings but a similar structure. We use standard nota- 
tion, Ck , where C denotes the minimal number of crossings 
in a projection onto a plane, and k distinguishes between 
knots with the same number of crossings [25] . 

To study the behaviour of knots in polymers under ten- 
sion we used a bead-spring model coupled to a coarse 
grained solvent [12l[26l[23, briefiy described here. The 
beads interact via the potential: 
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where the first term is an excluded volume interactions 
and the second term is a FENE spring potential to sim- 
ulate bonds. The parameters were chosen as: /3e = 1, 
a = 1, fc = 30 and Rq — 1.5. The equations of motion 
were integrated with a velocity Verlet algorithm. 

In this paper we used flexible polymers as we are mainly 
interested in generic, qualitative effects. For a quantitative 
comparison to experiments, it would be important to also 
include a bending potential that captures the semi-flexible 
nature of DNA, as has been done in previous simulations 
[2T|I24) . For example, in the elastic regime we study, semi- 
flexibility can lead to knots that extend fewer persistence 
lengths [20] than those in our simulations. 

The polymer was coupled to a mesoscopic solvent mod- 
elled with stochastic rotation dynamics (SRD) 28 . On 
average there are 5 solvent particles per cr^, and they 
provide a thermostat which conserves momentum and so 
preserves hydrodynamic interactions between monomers. 



The model thus describes a flexible polymer undergoing 
Brownian motion in a hydrodynamic solvent |261I27| . 
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Fig. 2: The two geometries used in the simulations. In (a) 
the polymer ends are connected using periodic boundary con- 
ditions, forming a ring polymer under tension. In (b), the 
polymer is linear and the two end beads are subjected to a 
constant force to provide the tension. Grey circles indicate the 
beads identified as the first and last beads of the knot by the 
knot tracking algorithm 

To simulate a polymer under tension, two different ge- 
ometries, illustrated in Fig. [21 were considered. In the 
geometry (a), a simulation box of 32a by 32a by 20a or 
21cr was used. The polymer was connected to itself across 
the periodic boundary. 

By varying the number of beads in the polymer, the 
tension can be changed. This approach offers the advan- 
tage of no free ends. Simulations may be run as long as 
necessary for good statistics without worrying about un- 
knotting or other end effects. In the geometry (b), a much 
larger simulation box of 32cr by 32a by 300cr was used and 
the polymer ends were not joined. Instead the first and 
last beads were subjected to a constant force. Geome- 
try (b) was mainly used to test the reliability of the more 
efficient geometry (a). 

Knots were introduced by hand. The first and last beads 
of a knot, shown in Fig. [21 can be identified by finding bond 
crossings of bead to bead vectors from the projection in 
the x-y plane (see Fig. [2]). The midpoint between the first 
and last beads was taken as the knot's position and the 
difference as its length. Distances were measured along 
the polymer contour, as in the experiments [20] . 

Fig. [3] shows that the length of the 3i knot decreases 
with tension and plateaus at higher tension, consistent 
with ref [21]. This behaviour holds for other knots too. 
We will primarily use a tension of / = 5 
regime since for a flexible polymer cr ~ /p. 
modeling DNA, where Ip sa 50nm, our force would map 
onto / ~ 0.4 pN, which is comparable to the experiments 
of Bao et al. [20], which were also in the elastic regime. 

In Table [H we compare the average length of different 
knots for polymers under a tension of / = 5 -^^^ . The knot 
length increases approximately linearly with the number 
of essential crossings, an effect we observe at other ten- 
sions as well. For geometry (a) it is harder to set tension 
explicitly. Instead either the number of polymer beads was 



, in the elastic 
If we were 



p-2 



Effect of topology on dynamics of knots in polymers under tension 



50| . . . . 1 

40- T 

r 30- T 

- % 

o 20- -'-f 

10 * 

°0 5 10 15 20 25 

Tension (kgT/o) 

Fig. 3: Average knot length for the 3i knot on a polymer under 
tension. Simulations are done with geometry (b). 

Table 1: Average knot length for linear polymers under a ten- 
sion of 5^^^ for geometries (a) and (b). 

Knot type Length in (b) Length in (a) 
geometry(cr) geometry(cr) 



3i 17 ±1 17 ±1 

4i 24 ± 2 24 ± 1 

5i 28 ±2 29 ±2 

6i 35 ±2 35 ±2 

7i 39 ± 2 39 ± 2 

8i 46 ± 3 45 ± 2 

9i 49 ± 3 50 ± 2 

lOi 56 ±3 57 ±2 

111 59 ±3 59 ±2 

12i 66 ±3 67 ±2 



changed, or the box length was varied between L — 20a 
and L = 21a (the SRD algorithm uses boxes of size a so 
constrains us to integer values). By comparing the knot 
length to the linear geometry, a similar effective tension 
could be simulated. We tested that the knot diffusion co- 
efficient only depends very weakly on tension in the range 
/ = 2 - , as expected for the elastic regime [20l[24] . 

We begin by studying the diffusive motion of the 
knots along the chain. In Fig. |4] we show the mean 
square displacement of the torus and even-twist knots as 
a function of time, measured in simulation units to — 
{a/2)y^mf /ksT, where m/ is the mass of a fluid parti- 
cle. At shorter times, we observe sub-diffusive behaviour 
with < >~ y/i, but this crosses over to Fickian dif- 
fusive motion, typically around t > lO^to. The resulting 
diffusion coefficients, extracted from the mean-square dis- 
placements at later times (typically t = 10^ — 2 x lO^to), 
are listed in Table [2] The diffusion coefficients clearly de- 
crease with increasing crossing number, but for a given 
complexity, torus knots have larger diffusion coefficients 
than even-twist knots. We observed a similar topology in- 
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Fig. 4: Mean squared displacement curves plotted on log-log 
scales. The upper lines are for the torus knots and have been 
displaced on the vertical axis by a factor of 10. The lower lines 
are the even-twist knots. Within each group knot complexity 
increases from top to bottom. The dashed lines have slopes of 
0.5 and 1. 



Table 2: Relaxation timescales and diffusion coefficients for 
knots in a polymer under tension. Timescales were found by 
an exponential fit to the autocorrelation decay between 0.5 and 
0.05, see Fig. [1 



Knot 
type 


Relaxation timescale 


Diffusion coefficient 

\ to ' 


3i 


0.020 ±0.001 


8.0 ±0.4 


4i 


0.029 ±0.001 


3.1 ±0.2 


5i 


2.1 ±0.1 


4.5 ±0.3 


6i 


1.2±0.1 


2.5 ±0.2 


7i 


2.5±0.1 


3.4 ±0.2 


8i 


1.4± 0.1 


2.0 ±0.1 


9i 


3.5±0.1 


2.3 ±0.2 


lOi 


1.9±0.1 


1.4 ±0.1 


111 


5.2±0.1 


2.1 ±0.1 


12i 


2.1±0.1 


1.3 ±0.1 



duced non-monotonic behaviour with crossing number in 
simulations of knotted polymer ejection from a capsid [T^ . 

The diffusion coefficient can be related to the friction 
coefficient through the Einstein relation ^ = ksT/D. We 
show results for friction vs. knot length in Fig. \5\ Both 
knot groups show the same linear dependence on length, 
with the same slope, but there is a constant shift between 
them. A very similar linear dependence of the friction 
on length was measured in experiments [20] as well as in 
simulations [24]. In both cases the only even- twist knot 
measured was 4i, and the friction was consistently higher 
than that expected from torus knots. We also plot a simple 
approximation to the friction, the average knot length, 
where ( is the friction coefficient of a monomer in the 
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Fig. 5: Friction coefficient against average knot length for torus 
knots (crosses) and even- twist knots (circles). The solid line is 
average knot length, where is the friction coefficient of a 
monomer in the SRD solvent (each bead contributes a length 
of cr.) 



SRD solvent. This simple approximation provides a good 
estimate of the slope of the friction vs. knot length curve 
in this tension regime, as also seen in ref. |24) . 

For miknotted polymers, neglecting the hydrodynamic 
interactions (Rouse approximation) leads to a linear de- 
pendence on the number of beads, but taking the hy- 
drodynamics into account (Zimm approximation) changes 
the scaling substantially [13 . The methods used in 
refs. |211[24] neglect hydrodynamic interactions, but the 
simulation technique we use here can reproduce the cor- 
rect Zimm scaling for linear polymers [1S1I17]- Includ- 
ing hydrodynamic interactions quantitatively changes the 
magnitude of the knot diffusion coefficients by about a 
factor of 1.5 but does not appear to alter their qualitative 
scaling with knot length. 

The difference in friction between torus and even-twist 
knots can be linked to an important topological distinc- 
tion. As knots diffuse along the chain, they must reptate 
through themselves [20]. From Fig. [1] one can appreciate 
that the polymer may pass relatively smoothly through 
a torus knot, always curving in the same sense. For an 
even-twist knot, however, the twist at one end introduces 
a sharp inversion of the direction in which the polymer is 
curving, making the polymer passage less smooth and in- 
creasing the effective friction. Since this inversion occurs 
only once in an even- twist knot, it may explain why the 
frictions are higher than those of the torus knots by an 
additive constant. 

Evidence that the main mode of knot diffusion in this 
regime is through reptation can be found in the mean- 
square displacement data in Fig. HI At short times the 
motion is sub-diffusive because as the polymer in the knot 
tries to move it tends to be bounced back by other sections 
of the knot. At longer times, when the polymer contour 
has relaxed within the constraints formed by the knot, we 



observe ordinary Fickian diffusion. This crossover corre- 
sponds to two regimes in the reptation model [13] where 
the mean squared displacement of the polymer along its 
own contour increases with the exponents we observe. 
Our results are also not inconsistent with a "sliding knot" 
picture such as in ref. 

The mapping of dynamics from coarse-grained simula- 
tion units to physical units is always subtle [29j . In the 
experiments of ref. [5D], D « 1.25fim^ / s for a 3i knot. 
Comparing to our simulated value for the same knot, 
and assuming Ip ^ a gives a mapping of to ~ 16/isec. 
Very similar values for to are found when using the other 
knots to do the mapping. This suggests that sub-diffusive 
< x'^ >^ \ft behaviour should persist up to a timescales 
of order ms, with the full crossover to Fickian diffusion oc- 
curing at time scales on the order of s. These time scales 
correspond qualitatively to the experiments of ref. |20j . 




Fig. 6: Decay of the knot length autocorrelation with time. 
Results for torus knots are plotted using solid lines and for 
even-twist knots using dashed lines. The inset shows 3i and 4i 
results on a different scale. 

We next investigate the fluctuations and dynamic modes 
of the diffusing polymer knot under tension. Orlandini 
et al [16] showed that the long timescale originally ob- 
served in the radius of gyration autocorrelation function 
in the "blob" regime can also be observed in the 
in the decay of the knot length autocorrelation function 

<U^g'/(0)>:</(0)></(0)> ' ^he^e ^(*) is the length of the knot 
at a given time. We investigated the same autocorrela- 
tion function for knots under tension, and plot the results 
in Fig. |6l The decay of the knot length autocorrelation 
function is approximately exponential, and the associated 
timescale (the inverse of the slope of the lines) of this expo- 
nential decay depends on the knot type, see Table ID The 
timescale increases with knot complexity and is also sig- 
nificantly longer for torus knots than for even-twist knots 
of similar complexity. There is a remarkably large jump 
(a factor of about 100) in the decay timescale between 
the most simple knots (3i and 4i) and the more complex 
knots (5i and larger). Preliminary investigations of poly- 
mer length, N, dependence showed that the timescales for 
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the more complex knots vary weakly with N, whereas the 
3i and 4i timescales are more sensitive to N. 
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Fig. 7: (a) Examples of knot modes (for 7i). Arrows indicate 
the direction of the tension force. The change between the left 
and right involves inversion of the spiral - the spiral-flipping 
mode. The change between top and bottom involves switching 
the ahgnment of the spiral with respect to the tension force - 
the breathing mode, (b) Origin of the faster breathing mode 
in even- twist knots. Left: in even- twist knots, the end of the 
spiral, the twist, may slide relatively freely along the strand 
that passes through it. Right: if the end of the spiral moves in 
a similar manner in a torus knot it collides with other loops. 

To identify the source of the slow decay of the knot 
length autocorrelation function, the knot configurations 
were more closely investigated. Two different modes were 
identified: a spiral-flipping mode and a breathing mode, 
both depicted in Fig.[71[a). 

The spiral-flipping mode involves switching between a 
state where the strand entering from the left predomi- 
nantly curves around the strand entering from the right 
and vice- versa. The breathing mode involves switching be- 
tween a state where the axis of the spiral is aligned with 
the tension force and one where the spiral forms a loop 
perpendicular to it. To characterise the knot state, two 
order parameters were used: one for each mode. 

For the spiral-flipping mode the sum of angles between 
successive bonds - in the plane defined by those bonds 
- multiplied by the displacement from the centre of the 
knot was taken over all of the beads in the knot. Depend- 
ing on whether the polymer is more curved towards one 
end of the knot or the other, this quantity is positive or 
negative. For the breathing mode the ratio between the 
maximum distance between any two beads in the knot in 
the direction of the tension force to the maximum distance 
perpendicular to the tension force was taken. When the 
spiral of the knot is aligned with the tension this quantity 
is larger. 

We calculated the autocorrelation functions, 

<0(t)0(O)>-<0(O)><0(O)> ^>,„^„ A :„ thp nrdpr na 

<0(o)0(o)>-<</.(o)><0(o)>' wnere (p is tne order pa- 
rameter value, for both order parameters. As an example. 
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Fig. 8: Autocorrelation functions for ffi for the knot length 
(solid line), the breathing mode order parameter (dashed line) 
and the spiral-flipping mode order parameter (dotted line). 
The inset shows results for 3i on a different scale. 

the autocorrelation functions for the 3i and 11 1 knots 
are depicted in Fig. [Sj The breathing autocorrelation 
function has a very similar decay to the knot length 
autocorrelation function, whereas the spiral flipping au- 
tocorrelation function decays more rapidly. Nevertheless, 
both modes show a large jump in time scale between the 
3i and 4i and the more complex knots. 

These results suggest that the breathing mode is re- 
sponsible for the slow decay of the knot length autocorre- 
lation. The average knot lengths for high and low values 
of the breathing mode order parameter support this. The 
extended configuration has a larger knot length; for the 
most complex knots a difference of about 3cr is seen. 

The coupling to the breathing mode offers an explana- 
tion for the large jump in the timescale between 3i and 
4i and the more complex knots. The spiral of the most 
simple knots consists of only one loop and so shifting from 
a configuration which is extended in the tension direction 
only involves changing its shape. In contrast, for more 
complex knots, it involves co-ordinated movement of mul- 
tiple loops. It also offers an explanation for the differ- 
ence between torus and even- twist knots. As depicted in 
Fig. El^b), an even- twist knot may deform from being ex- 
tended by sliding the twist along the strand that passes 
through it. In a torus knot the equivalent movement is 
blocked by the loops of the spiral so breathing is slower. 

Note that there is no big jump in the diffusion coef- 
ficients between 3i and 4i and the more complex knots. 
While the breathing and spiral-fiipping modes may affect 
the crossover from sub-diffusive to diffusive behaviour (in 
fact the sub-diffusive range is considerably shorter for the 
3i and 4i knots) they do not appear to determine the long- 
time diffusive dynamics. Further evidence follows from the 
fact that the torus knots have slower relaxation but faster 
diffusion than the even- twist knots. 

It is instructive to compare our results for the longest 
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relaxation timescale to ring polymers without tension in 
ref. [TS]. By contrast with our results, the long timescale 
was shown to decrease with increasing knot complex- 
ity and there was no significant difference between knot 
groups. This suggests that the dominant modes determin- 
ing relaxation in the "blob" regime are not related to the 
breathing modes we observe for polymers under tension. 

By contrast, simulations in ref. [18] that studied the re- 
laxation of initially knotted polymers that are cut (with- 
out tension) exhibited a relaxation timescale that in- 
creased with knot complexity, and was longer for torus 
knots than for even-twist knots. The similarity is perhaps 
not surprising. For a cut knot to become unentangled the 
polymer must move along the contour formed by the erst- 
while knot, a process analogous to the self-reptation [50] 
of knots under tension. 

In summary: we exploited a new boundary condition 
that minimizes end-effects to study the effect of topology 
on the dynamics of knots in polymers under tension in the 
"elastic regime" where the force is strong enough to affect 
the behaviour of the knot, but not so strong that details 
of the inter-molecular potentials dominate the behaviour. 

The observed motion of the knot, including a sub- 
diffusive regime of the mean-square displacement at short- 
times, is consistent with a self-reptation mechanism ^0\. 
To first order, the friction coefficient scales linearly with 
knot length. Topological differences also play a role. For 
a similar knot length or complexity, torus knots have a 
smaller friction (and so a larger diffusion coefficient) than 
even-twist knots. 

Fluctuations and dynamic modes are also affected by 
knot topology. The longest relaxation mode of the length- 
length autocorrelation function is dominated by a breath- 
ing mode. This topological timescale increases with knot 
complexity and is significantly longer for the torus knots 
than for even-twist knots of similar complexity. These dif- 
ferences can be rationalized by examining knot topology. 

In this paper we have used a simple flexible polymer 
model. The main results summarized should carry over 
to other polymer models. Nevertheless, it would be in- 
teresting, for example, to simulate semi-flexible polymers 
that provide a better DNA model. In particular, this may 
give a more accurate estimate of physical timescales, and 
aid experiments to measure them. In addition, it would 
be important to study the role of tension on DNA knots 
in vivo [23]. What are the magnitudes of the biological 
forces, and what is the role of macromolecular crowding 
on knot diffusion and relaxation rates? Finally, the advent 
of nanotechnology has brought questions about the diffu- 
sion and dynamics modes of knotted polymers in strong 
confinement to the fore. It would be interesting to ex- 
tend work that looks at the interplay of confinement and 
tension [301131] . 

^ * * 

We thank Olivier Pierre-Louis and Enzo Orlandini. 
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